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ABSTRACT 

Re-observations with the H.E.S.S. telescope array of the very-high-energy (VHE) source HESS J1018-589 A coincident with the 
Fermi-LNI y-ray binary 1EGLJ1018.6-5856 have resulted in a source detection significance of more than 9a, and the detection 
of variability ix^jv of 238.3/155) in the emitted y-ray flux. This variability confirms the association of HESS J1018-589 A with the 
high-energy y-ray binary detected by Fermi-LNI, and also confirms the point-like source as a new very-high-energy binary system. 
The spectrum of HESS J1018-589 A is best fit with a power-law function with photon index T = 2.20 ± O.ldjtat ± 0.2sys. Emission is 
detected up to ~20TeV. The mean differential flux level is (2.9 ± 0.4) x 10“'^ TeV^'cm^^s^' at 1 TeV, equivalent to ~1% of the flux 
from the Crab Nebula at the same energy. Variability is clearly detected in the night-by-night lightcurve. When folded on the orbital 
period of 16.58 days, the rebinned lightcurve peaks in phase with the observed X-ray and high-energy phaseograms. The fit of the 
H.E.S.S. phaseogram to a constant flux provides evidence of periodicity at the level of N„. > 3a. The shape of the VHE phaseogram 
and measured spectrum suggest a low inclination, low eccentricity system with a modest impact from VHE y-ray absorption due to pair 
production (t < 1 at 300 GeV). 

Key words, gamma rays: stars; X-rays: binaries; stars: individual: IFGLJIOI 8.6-5856; acceleration of particles; radiation mechanisms: 
non-thermal 
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1. Introduction 


The region around the supernova remnant (SNR) SNRG284.3- 
1.8 (Milne et al. 19891 shows two clearly distinct regions 


of very-high-energy (VHE; E>100GeV) gamma-ray emis¬ 
sion ( |Abramowski et aT]|2012| l; an extended emission named 
HESS J1018-589 B likely associated with the pulsar wind nebula 
(PWN) powered by the bright pulsar PSR J1016-5857 (|Camilo| 
jet al.|2001|[2004| , and the point-like source HESS J1018-589 A. 
The latter is positionally coincident with IFGLJIOI8.6-5856, a 
point-like high-energy gamma-ray (HE; 100MeV<E<100GeV) 
variable source detected by the Fermi Large Area Telescope (EAT) 
dAbdo et al.|20T0l i. 

The y-ray binary 1FGLJ1018.6-5856 was detected in a blind 
search for periodic sources in the Fenni-LAJ survey of the Galac¬ 
tic Plane through the modulation of its HE y-ray flux ( |Acker-| 
mann et al. 2012|i. Optical observations show that the non-thermal 


source is positionally coincident with a massive star of spectral 
type 06V((f)). The radio and X-ray flux from the source are 
modulated with the same period of 16.58+0.02 days, interpreted 
as the binary orbital period (jPavlov et al.|2011[ Li et al.|201 


[Abramowski et al.|2012[|An et al.|2013| ). 

The spectrum of the periodic source in the Fermi-hPfF do¬ 
main exhibits a break at ~ 1 GeV with best-fit values of rHE(0.1-1 
GeV)=2.00+0.04 and rHE(l-10 GeV)=3.09+0.06 and an inte¬ 
gral energy flux above 100MeV of (2.8+0.1)xl0“'‘’ergcm“^ s“'. 
The HE y-ray spectral shape evolves with orbital phase, with a 
decrease in spectral curvature at flux minimum of the emission 
(between phases 0.2 and 0.6) and a hardening of the spectrum at 
flux maximum. 

The best-fit position reported in the 
(| Abramowski et 3^1 |2012|) 

nn- - 


per 


previous pa- 
for HESSJ1018-589A, 


a=10"18”59.3'^±2.4^(3( and (5=-58°56'10"±36"^t (J2000), is 


compatible with the 95% confidence contour of IFGL J1018.6- 
5856. The VHE emission is well-described by a power-law 
function with a spectral index of r=2.7±0.5stat±0.2sys, similar 
to the one describing the VHE emission of the larger region 
HESS J1018-589 B (r=2.9±0.4stat±0.2sys). No variability was 
found in the H.E.S.S. data set, although the contamination of 
the nearby source and the uneven sampling of the observations 
prevented a Arm conclusion at the time. 

Here, a deeper study of HESS J1018-589 A to assess its asso¬ 
ciation with the y-ray binary is presented. In Section]^ the data 
sample and results are described. In Section the features of 
HESS J1018-589 A are discussed in light of the multi-wavelength 
observations available, and conclusions are drawn in Section]^ 


2. Data Analysis and Results 

The H.E.S.S. telescope array is a system of five VHE y-ray imag¬ 
ing atmospheric Cherenkov telescopes (lACTs) located in the 
Khomas Highland of Namibia (23° 16'18" S, 16°30'00" E). The 
fifth telescope was added to the array in summer 2012 during the 
H.E.S.S. phase-II upgrade, increasing the energy coverage and 
boosting the system sensitivity. The nominal sensitivity of the 
H.E.S.S. phase-I array (excluding the large telescope) reached 
in 25 hours is ~2.0xl0^'^ ph cm^^s ' (equivalent to 1% of the 
Crab Nebula flux above 1 TeV) for a point-like source detected 
at a significance of 5cr at zenith. The stereoscopic approach re¬ 
sults in a positional reconstruction uncertainty of ~6' per event. 
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Fig. 1. The SED of HESS J1018-589 A / IFGL J1018.6-5856 is shown 
in black (filled squares and circles for the LAT and H.E.S.S. detection 
respectively). Eor comparison, the SEDs of LS 5039 during superior 
(SUPC) and inferior conjunction (INFC) are also included (blue points 
from [Hadasch et al.|20f^ [Aharonian et al.|2005a) 


good energy resolution (15% on average) and an efficient back¬ 
ground rejection ( [Aharonian et al. 2006) 1. H.E.S.S.-I observed 
the region towards the Carina arm from 2004 to 2009. The data 
set presented in [Abramowski et’ak ( 2012| l was increased in the 
subsequent years from 40 hours to 63.3 hours effective time with 
dedicated observations with the H.E.S.S.-I array to cover the 
orbital phases in which the emission of IFGL J1018.6-5856 ob¬ 
served in X-rays and HE increases. The zenith angle at which 
the source was observed ranged from 35° to 55° resulting in a 
mean energy threshold of 0.35 TeV. These new observations were 
performed in wobble-mode, during which the telescopes were 
pointed offset (0.7°) from the nominal source location to allow si¬ 
multaneous background estimation. The data were analyzed using 
an improved analysis technique (multivariate analysis, reaching 
~0.7% of the Crab Nebula flux at 1 TeV, 5 <t; [Becherini et all 
[2011[ | and cross-checked with the Hillas second moment event re- 
construction method ([Aharonian et al.|2006) l and Model Analysis 
( de Naurois & Rolland|2009[ l, including independent calibration 
of pixel amplitudes and identification of problematic or dead pix¬ 
els in the lACTs cameras. The spectra and light curves shown 
here are derived for a cut of 80 photoelectrons in the intensity of 
the recorded images. 

The new analysis of HESS J1018-589 A, using the larger 
data set, confirms the point-like VHE y-ray emission reported in 
Abramowski et al.|([2012)i. The y-ray signal is detected with a sta¬ 


tistical significance of 9.3cr pre-trials (derived using an oversam¬ 
pling radius of 0.10° and corresponding to more than 1.5cr post¬ 
trials), centered at Q'=10M8”58*±5^,^t and 5=-58°56'43"±30",, 
(J2000). The best-fit position is estimated by means of a maxi¬ 
mum likelihood fit of the exposure-corrected uncorrelated excess 
image. This position is compatible with the position derived in 
[Abramowski et al.[ ( 2012[ i, but the presence of the nearby extended 
source HESS J1018-589 B precludes from an improvement in the 
position uncertainty even with the additional observation time. 
The fitted extension is compatible with the H.E.S.S. point spread 
function (PSF, estimated to have a mean 68% containment radius 
of ~0.1°). The obtained position is used to derive the spectrum 
of the point-like source, integrating in a circle of 0.1° around 
it and using a forward-folding maximum likelihood fit ([Piron 


et al.|200r i. The photon spectrum is well-described with a power- 


law function with index r=2.20+0.14stat±0.2sys (Fig.[TJl and the 
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Fig. 2. Lightcurve of the integral flux above 0.35 TeV in a 0.1° region 
centered on HESS JIO18-589 A binned by observation run, correspond¬ 
ing to approximately 30 minutes of observation time per bin. The dashed 
horizontal line shows the mean integral flux. 
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flux normalisation is Nq — (2.9 + 0.4stat) x 10“^^ TeV“ 
at 1 TeV. The systematic error on the normalisation constant 
No is estimated to be 20% ( jAharonian et al.|[20061 l. The better 
statistics allow for a better determination of the spectral fea¬ 
tures of the point-like source compared with the one presented 
1 Abramowski et al.|(|2012|l, including a clearer separation from 


HESS J1018-589 B. The nearby source introduces a maximum 
of 30% contamination on HESS J1018-589 A, although above 
1 TeV, thanks to the better PSF, less than 10% contamination was 
calculated from a simultaneous fit of the two sources. 


The light curve of the source above 0.35 TeV, binned by ob¬ 
servation run (approximately 30 minutes of observation time), is 
shown in Eig.|^ The best-fit mean flux level above 0.35 TeV is 
marked with a dashed gray line. The lightcurve displays clear vari¬ 
ability, with ax^lv of 238.3/155 (corresponding to 4.3cr) using a 
likelihood ratio test with a constant flux as null hypothesis. 

To investigate the periodicity of the source, the data were 
folded with the 16.58 day period found in the HE y-ray ob¬ 
servations (Eig. 1^ top panel) using the reference time of 
Tmax=55403.3 MJD as phase 0 ( |Ackermann et al!]|2012 i in a 
single trial. The number of bins in the phaseogram was selected 
to obtain a significance of at least Icr in each phase bin. For 
comparison, the same phaseogram is also shown for HESS J1023- 
589, a nearby bright y-ray source expected to be constant. The 
flux variation along the orbit shows a similar behaviour when 
comparing it with the Fermi-LA^T flux integrated between 1 and 
10 GeV (Fig.[^ middle panel). An increase of the flux towards 
phase 0 is observed, with a^^/v of 22.7/7 (3. Icr) when fitting the 
histogram to a constant flux, providing evidence of periodicity at 
the a priori selected period. Unfortunately, the uneven sampling 
and large timespan of the observations did not allow for an inde¬ 
pendent determination of the periodicity from the VHE y-ray data 
using a Lomb-Scargle test ( |Scargle|1982| l, since the equivalent 
frequency is ~8 times larger than the sample Nyquist frequency. 
Finally spectral modulation was examined by deriving the photon 
spectrum for observations in the 0.2 to 0.6 phase range (motivated 
by the Fermi-LAH observations) and comparing it with the one 
derived at the maximum of the emission in the complementary 
phase range. No spectral modulation was found within the photon 
index errors (AT = 0.36 + 0.43) although it should be noted that 
the data statistics in the 0.2 to 0.6 phase range are insufficient (3cr 
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Orbital Phase 

Fig. 3. VHE, HE, and X-ray fluxes of lEGL J1018.6-5856 folded 
with the orbital period of P=16.58d. Two orbits are shown for clar¬ 
ity. Top: VHE integral flux above 0.35 TeV measure by H.E.S.S. (red 
circles). Eor comparison, a scaled lightcurve from the nearby bright 
source HESS J1023-589 is shown in gray. Middle top and middle bottom: 
Fermi-hAI lightcurve between 1 and 10 GeV (solid blue squares) and 
between 0.1 and 1 GeV (open blue squares; [Ackermann et al.|2012[ l. 
Bottom: X-ray 0.3-10 keV count rate lightcurve from 67 Swi^-XRT 
observations in 2011 (green), 2012 (blue), and 2013 (red). 


detection) to firmly conclude a lack of variation in the spectrum 
at different orbital phases. 

In order to compare the VHE orbital modulation with the 
behaviour of the source at X-ray energies, 67 Swift-XRT obser¬ 
vations of 1FGLJ1018.6-5856, performed between 2011 and 
2013 and with a median observation time of 2.2ksec, were anal¬ 
ysed. Early subsets of these observations were presented previ¬ 
ously by [Ackermann et ^ ( |2012[ ) and jAn et al.| ( |2013| l. Cleaned 
event files were obtained using xrtpipeline from HEAasoft 
v6.15.1. For each observation, source count rates were extracted 
from a 1 arcmin circular region around the nominal position of 
1FGLJ1018.6-5856, and background count rates extracted from 
a nearby region of the same size devoid of sources. The resulting 
count rate lightcurve, folded with the orbital period, is shown in 
the bottom panel of Figure]^ The phaseogram displays a sharp 
peak around phase 0, matching the location of the maximum in 
the VHE and HE phaseograms. There is an additional sinusoidal 
component with a maximum around phase 0.3 and with lower 
amplitude than the sharp peak at phase 0. 
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3. Discussion 


The flux variability and periodical behaviour of HESS J1018- 
589 A suggest the identification of the VHE source with the 
y-ray binary 1FGLJ1018.6-5856. It therefore becomes the fifth 
binary system, along with LS 5039 (jAharonian et al.||2005a|), 
LSI+61 303 dAlbert et al.||20061l, PSRB1259-63 (|Aliaronian| 
let al.|2005b) and HESS J0632+057 dAcciari et al.|2009| l detected 
at VHE during multiple orbits, in addition to the hint of a flaring 
episode from the X-ray binary Cygnus X-1 dAlbert et al.|2007 1. 
When folded with the modulation period found at other wave¬ 
lengths, the rebinned VHE lightcurve shows a modulation, signif¬ 
icant at the 3. Icr level, in phase with the HE y-ray lightcurve. The 
phaseogram (Fig.[^ shows a similar behaviour (within the limited 
statistics) to the high-energy lightcurve of IFGL J1018.6-5856 
showing a flux increasing simultaneously to the one occurring in 
the HE and X-ray counterpart. 

Despite the different orbital behaviour at different wave¬ 
lengths, the stars in the y-ray binaries IFGL JlOl8.6-5856 and 
LS 5039 are thought to be very similar, with spectral types of 
06V((f)) and 06.5V((f)), respectively dAckermann et al.|2012[ 
jClark et al.||200l| . Unfortunately, the orbital parameters of 
HESS J1018-589 A are not yet known and only limited conclu¬ 
sions can be drawn on the relation between the compact object and 
the massive star. Both binary systems are composed of an almost 
identical massive star and a compact object orbiting around it on a 
timescale of days. The four-times larger period of IFGL J1018.6- 
5856 implies a factor ~2.5 larger semi-major axis than in LS 5039, 
and the low amplitude of the flux modulation observed by Fermi- 
LAT, of the order of 25%, can be interpreted as a sign of a low- 
eccentricity orbit. Although the spectral index variability at HE 
y-ray is at odds with anisotropic IC being the only source of flux 
variability, such low modulation amplitude would be difficult to 
realize under the widely changing conditions of an eccentric orbit. 
The behaviour of HESS JlOl8-589 A at different orbital phases 
is mimicked in X-rays, HE, and VHE, showing in all cases a max¬ 
imum flux near phase 0. There is a second sinusoidal component 
that peaks at phase 0.3 and appears in radio ( |Ackermann et al.] 
2012|l and X-rays (Fig.|^ bottom panel), as well as a hint in the 


0.1 to 1 GeV Fermi-LPCf lightcurve peaking at phase 0.5 (Fig.|^ 
middle-bottom panel), but it is not observed (with the current 
statistics) at higher energies. However, in LS 5039 the VHE flux 
is correlated to the X-ray flux but anti-correlated to the HE flux 


(|Aharonian et al.|2005at Takahashi et al.|200^|Hoffmann et al!] 


2009[ Hadasch et al.|2012^ . The HE spectral energy distributions 

(SEDs) of the two binary systems are remarkably similar in shape 
and flux (see Fig. [^l, although it should be noted that the systems 
are believed to be located at different distances: whereas LS 5039 
is ~2.5 kpc away, IFGL J1018.6-5856 is believed to be located 
at 5+2 kpc, derived from the interstellar absorption lines of the 
companion ( Ackermann et al.|2012[ ). At VHE, LS 5039 shows a 
clear spectral modulation at different orbital phases, with a mean 
luminosity between 1 and 10 TeV of ~ 10 ^^(£/ 2 . 5 kpc)^ergs^\ sim¬ 
ilar to the one found in HESS JlOl 8-589 A in the same energy 
range (9.9 x 10 ^^((i 5 kpc)^ergs"'). However, the ratios between 
the fluxes measured at HE and VHE of the two binary systems 
differ substantially (see Figs.[T]an d|3): whereas for LS 5039 the 
ratio between the fluxes at 1 GeV and 1 TeV varies between ~15 
and 40 in superior and inferior conjunction, respectively, for 
1FGLJ1018.6-5856 and HESSJ1018-589 A aratio of-160 is 
found, with the TeV flux strongly reduced with respect to the 
GeV flux when compared with LS 5039. Similar to other binaries, 
the spectrum of IFGL J1018.6-5856 measured at HE does not 
extrapolate to VHE. 
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Fig. 4. Results of the fit of a pair-production-absorbed power-law to the 
HE and VHE spectra of fFGLJfOl8.6-5856, with solid (dashed) lines 
indicating the absorbed (intrinsic) spectra. The emitter is assumed to be 
in the plane of the sky with respect to the star (i.e., on the plane of an 
orbit with i = 0°). The model with the maximum absorption compatible 
with the VHE spectral data is shown in red, the best-fit model with 
(rsoo) = 1 in green, and the best-fit model with no absorption in blue. 
Observational points are as in Fig.[T] 


Regardless of the nature of the emission process responsible 
for the TeV emission, the strong stellar photon field in the envi¬ 
ronment of the binary system unavoidably leads to the absorption 
of y-rays above -50 GeV through pair pro duction (|Moskalenko &| 
|Karakula|1994|[B6ttcher & Dermer|2005HDubus|2006| ). Assum- 

ing that the VHE emission is due to anisotropic inverse-Compton 
(IC) scattering of a leptonic population at a similar location as the 
HE one, the identical phasing of HE and VHE would imply a low 
inclination of the orbit with respect to the observing direction. 
Furthermore, the sinusoidal modulation implies that the orbit 
cannot be highly eccentric (or the density of photons would lead 
to strong variation). Under this scenario, IC anisotropic emission 
would be most efficient when the emitter is behind the star with 
respect to the emitter, i.e., at superior conjunction, and therefore 
the flux maximum at phase 0 can be tentatively associated with 
this orbital configuration. 

In order to illustrate the effects of absorption on the observed 
VHE spectrum, in the following calculations a circular orbit on 
the plane of the sky (/ = 0 °) is assumed to exemplify the low, but 
likely non-zero, inclination of the orbit, and take stellar parame¬ 
ters as in LS 5039. A fit of a pair-production-absorbed power-law 
function (of shape cx exp(- (t^)), where {te) is the energy- 
dependent orbit-averaged optical depth and a the intrinsic spectral 
index) to the measured VHE spectrum indicates that the maxi¬ 
mum optical depth at 300 GeV compatible with the VHE data 
(at 68 % CL) is (T 300 ) - 3.4, with an intrinsic index of a ^ 3.1. 
Figure]^ shows how the energy dependence of pair-production 
absorption results in a power-law-like spectrum between 500 GeV 
and 10 TeV even for high optical depths, as long as the intrinsic 
spectrum is steep enough. This means that for an optical depth 
of (T 300 ) - 3.4, the steep spectral index required to fit the VHE 
data would result in a large HE emission below 100 GeV, where 
pair-production absorption is no longer significant, up to a factor 
100 brighter than the flux observed by Fenni-LPfF between 10 
GeV and 100 GeV. Therefore, either the intrinsic emission from 
the VHE component has a sharp spectral break between 80 and 
200 GeV, or the VHE instrinsic spectrum must be significantly 
harder than a ^ 3. Considering the latter, and taking the Fermi- 
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LAT flux between 10 GeV and 100 GeV as an upper limit to the 
emission of the VHE component at these energies, the optical 
depth should be lower than 1, as illustrated by the green model in 
Fig. a For an orbital inclination of i = 0°, an optical depth lower 
than unity hints towards an emitter located farther away from the 
star than the compact object (at a distance of at least ~ 3x10'^ cm 
from the compact object). At higher orbital inclinations, the limit 
placed on the orbit-averaged optical depth can not be directly 
related to the location of the emitter, given that the optical depth 
would vary signiflcantly along the orbit. However, the correlation 
between HE and VHE emission and the sharpness of the peak of 
VHE emission at superior conjunction indicate that the optical 
depth at this position must be low enough to not have a signiflcant 
effect on the observed flux modulation, therefore excluding an 
emitter close to the compact object for high orbital inclinations. 


Several mechanisms have been proposed to explain VHE 
variability and periodic modulation via either IC processes or 
pion production of high energy protons with the companion wind 


(|KirketaI.|I999 

|Bosch-Ramon et aL|2006||Dermer & Bottcher| 

2006[ Bednarek 

|2006| Dubus 2006 |KhanguIyan et al.| 

2008 


[Sierpowska-Bartosik & Torres|2008| r. In a leptonic scenario, the 
maximum energy of the H.E.S.S. measured spectrum can be used 
to derive further constraints on the location, magnetic field and 
acceleration efficiency of the VHE emitter in HESS J1018-589 A. 
Given the energy of the stellar photons, IC scattering will take 
place in the deep Klein-Nishina (KN) regime, in which all of 
the electron energy is transferred to the scattered photons. In this 
scenario, the maximum energy detected (up to ~20 TeV) would 
require the presence of 20 TeV electrons in the VHE emitter, 
which in turn requires that they are accelerated faster than their 
radiative energy loss timescale. The acceleration timescale can 
be expressed as: 

Gcc “ ^acc^hl^ ^ 0.1 ^accExevBQ S, (1) 


where is the Larmor radius of the electron, Ejev is the elec¬ 
tron energy in TeV units, Bq is the strength of the magnetic held 
in Gauss, and ?7acc > I is a parameter that characterizes the ef¬ 
ficiency of the acceleration (in general » I, and only for 
extreme accelerators does rj^cc approach I, i.e. the Bohm limit). 
The balance between face and the cooling tir ne of electrons in 
the KN regime, given by fKN 


lO^c/jjE^gy s (|KhanguIyan et al. 


20081 , where di^ is the distance to the optical star in units of 


10'-’ cm, implies E„ 


(lO'^Bcflacc'^rlTeV. For IC domi¬ 


nant losses, and considering the maximum energy in the VHE 
spectrum, a minimum B > 2.5 x 10 "*77acc^fi3 G can be derived. 
Furthermore, if non-radiative (adiabatic) energy losses are neg- 
liglible, electron energy losses in the energy band relevant to 
VHE emission would be dominated by the interplay between 
IC losses, which in the KN regime decrease with energy, and 
synch rotron losses, which increase with energy ( |Moderski et al.| 
2005 I. For a power-law injection spectrum with canonical 
Pin] = 2, this results in a hardening (pe ~ 1-3) of the spectrum 
of the underlying steady-state particle population up to the en¬ 
ergy for which IC and synchrotron losses are balanced, and a 
softening (pe ~ 3) for higher energies (see, e.g., |Moderski et ^ 
2005 [ Dubus et al.|2008|l. The energy of the cooling break, Ebreak, 


can be found from the balance of IC and synchrotron cooling 
timescales fKN = Isyn, which, taking fj, 


TOOE^'y^G 


s, re¬ 


suits in Ebreak ~ 0.58 (Eg'^i 3) * TeV. The relatively hard VHE 
spectrum detected from 1FGLJ1018.6-5856 requires an evolved 
particle distribution with pe < 2, indicating that Ebreak should be 
higher, or of the order of, the electron energies sampled by the 
TeV spectrum. Considering Ebreak ^ 10 TeV, the magnetic field 
strength is constrained by the VHE spectrum to 5 < O.lc/jj* G. 


These constraints depend strongly on the location, acceler¬ 
ation efficiency and magnetic field of the emission region. An 
extended discussion of these relationships for a VHE emitter in 
a binary system can be found in [Khangulyan et al.| ( |2008) l. For 
the case of 1FGLJ1018.6-5856, the extension of a hard VHE 
spectrum up to 20 TeV indicates that acceleration/emission re¬ 
gions close to the compact object require an extremely efficient 
acceleration process, with rj^cc ^ 50 and magnetic held strengths 
between 0.001 and 0.1 G. If the emitter is located farther from 
the star, the constraint on the acceleration efficiency is relaxed, 
but the upper limit on the magnetic field is reduced to 0.03 G at 
d - 3 X 10'^ cm. Regardless of the location of the emitter, the re¬ 
quirement that the magnetic held strength is below 0.1 G indicates 
that, in this scenario, the particle population responsible for the 
VHE emission would have a maximum 2-10 keV X-ray ffux of 
1.2 X 10 '"^ erg cm“^ s', more than an order of magnitude lower 
than its detected X-ray ffux of (6.5 + 0.7) x 10 ergem^^ s ' 
(|Abramowski et al.|2012|), a similar situation to that found for 
LS 5039 dZabalza et al.|2013] l. 


4. Conclusions 

The new observations of HESS J1018-589 A with the 
H.E.S.S. telescope array have increased the significance of the 
detection up to ~9cr, allowing the firm identification of a new 
VHE binary system through the measurement of its variable emis¬ 
sion at a significance level of 4.3cr. Folding the measured flux 
on a 16.58-day orbit results in a phaseogram similar to the one 
observed at HE, with a wide peak around phase 0. The result 
of fitting the phaseogram to a constant flux indicates evidence 
of periodic flux at the 3.1cr level. The phase-averaged photon 
spectrum extends up to ~20 TeV, posing constraining limits on 
the magnetic held (0.001<B<0.1 G). Likewise, the spectral shape 
above 0.350 TeV limits the y-y absorption and optical depth at 
300 GeV to T(300GeV) < 1, which will be helpful in constrain¬ 
ing the location of the emitter once the orbital parameters are 
known. 

Deeper observations with H.E.S.S. II will improve the statis¬ 
tics at VHE and will provide a measurement of the spectrum 
below E < 100 GeV, allowing the investigation of key properties 
of the binary system such as spectral variation within the orbit 
or the spectral shape at low energies. Finally, the investigation 
of the orbital parameters through radio and optical observations 
is crucial to the understanding of the VHE emission mechanism 
in combination with the periodicity and variability observed at 
lower energies. 
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